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We study a square-lattice three-state Potts antiferromagnet with a staggered polarization field at 
finite temperature. Numerically treating the transfer matrices, we determine two phase boundaries 
separating the model-parameter space into three parts. We confirm that one of them belongs to the 
ferromagnetic three-state Potts criticality, which is in accord with a recent prediction, and another to 
the Ising type; these are both corresponding to the massless renormalization-group flows stemming 
from the Gaussian fixed points. We also discuss a field theory to describe the latter Ising transition. 
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It is widely recognized that strong frustrations can 
provide a ground state with an extensive entropy and 
prevent long-range orderings of systems. The simplest 
example may be the triangular-lattice antiferromagnetic 
(AF) Ising model, whose ground state is critical and ex- 
hibits power-law decays of correlation functions P, 
Although, unlike in the ferromagnetic (F) cases, the 
models of this type may depend on their details (e.g., 
lattice structures), the field theories used to describe 
the ground state and lower-energy excitations have at- 
tracted much attention. Further, recent interests in 
this area are rather focused on mutual relations of fixed 
points (e.g., crossovers of criticalities) embedded in the 
renormalization-group (RG) flows 0, so the understand- 
ings of their universal properties are quite important . 

The AF three-state {q = 3) Potts model on the square 
lattice A exhibits the same properties as those systems; 
its Hamiltonian is described by using the ternary vari- 
ables cFj = 0, 1, 2 {i e A) as 



(J>0), 



(1) 



where the sum runs over all nearest-neighboriiig (NN) 
pairs. While its ground-state P, IE H , H,M ^-^i^ 
finite temperature properties ll3| have been inten- 
sively investigated, a possibility of the crossover behav- 
ior from the AF to the F three-state Potts criticality has 
been recently proposed on the basis of a continuous field 
theory In close relation to the classical spin sys- 

tem, the realization of the field theory was found in the 
ID quantum spin system (i.e., a frustrated XXZ Heisen- 
bcrg chain in magnetic fields), and it can show the F 
three-state Potts criticality in its ground state In 
this Letter, we quantitatively investigate the possibility 
of crossover by introducing a square-lattice three-state 
Potts model defined below, which is also relevant to ID 
quantum systems. Before presenting the formulation of 
our investigation, we shall briefly refer to the relating 
research so far. 

Since the ground state of Hq is equivalent to the six- 
vertex model on the ice point, it shows the Gaussian 



criticality with the conformal anomaly number c = 1 
0, 0, 0, 13, 13 ■ This sort of equivalence was used to de- 
termine the properties of lower-energy excitations: den 
Nijs, Nightingale, and Schick found that both the uni- 
form Sj = g27r2drj/3 Q^^^ |.]^g staggered magnetizations 
Sj ~ {—lysj are relevant with scaling dimensions = | 
and xs = ^, respectively Here, (—1)-' = ±1 for j in 
the even (odd) sublattice A± . It was pointed out that the 
scaling dimensions of relevant scalar operators can take 
three values, and other than above two, the staggered 
polarization Pj = (— 1)"' X^U^'^tjj.o-fc — 1) takes = | 
[the sum is over k next-nearest neighboring (NNN) with 

j] i- 

Numerical studies have been also performed llll | : 
Salas and Sokal using the Monte Carlo (MC) simulations 
at zero temperature succeeded in confirming the lowest 
two scaling dimensions, but the estimation for the third 
one exhibits a deviation due to its larger energy scale 
The scaling dimensions can be also obtained by 
the transfer matrix technique; the lowest two were ac- 
curately obtained by de Queiroz, but the third could 
not be found 0|. More profound understanding at fi- 
nite temperature was brought about by Cardy, Jackob- 
scn, and Sokal ■ They pointed out that there are two 
types of excitations controlled by the thermal scaling field 
u = , i.e., the relevant one with = |, and the 

marginal one, and that both of these are necessary to 
explain the exotic corrections to scaling observed in the 
MC data 0|. Nevertheless, as its main role, the energy 
operator = X]fc ^crj,<yk (the sum is over k NN with j) 
brings about the second-order phase transition with the 
divergent correlation length ^ oc ^"^/(^-ajt) [^[f^. 

Recently, Delfino argued the "bosonization" descrip- 
tions of the transition and above-mentioned excitations 
. He also provided their discrete symmetry proper- 
ties, and then concluded that the criticality in the ground 
state can exhibit the crossover behavior to the three-state 
Potts criticality with c = | as a resultant of the compet- 
ing relevant perturbations Cj and Pj, and its description 
is given by the self-dual sine-Gordon model with the di- 
mcnsionlcss coupling 0^ = 67r (see also Ref. [T^). 
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Based on these developments, especially for the aim 
of quantitative understandings on the competitions be- 
tween relevant perturbations and the resultant crossovers 
of criticalities, we shall introduce the following model, 
namely, the square-lattice AF three-state Potts model 
with the staggered polarization field defined by the re- 
duced Hamiltonian H{K, V) = H/k-eT with 

n{K, V) = K ^ 5,^,,, -VY,{-IY6.^,.,. (2) 

The second sum runs over all NNN pairs [j, k] , so that 
the sublattice symmetry is explicitly broken for nonzero 
V , but the S3 symmetry associated with the global per- 
mutations of the ternary variables is preserved. Other 
than {K, V) = (00, 0), there are two special points in this 
model: When = 0, the system decouples into those 
defined on A±, and thus (0,^3'''') and (0,oo) correspond 
to the exact transition point of the F three-state Potts 
model on A+ and to the ground state of the AF three- 
state Potts model on A_, respectively \VI^^ := ln(f + ^)] 
[lij . Further, when X > 0, we expect that the disordered 
phase at y ^ 1 will change into a phase with the F or- 
dering on A-|_ and with the AF ordering on A_ satisfying 
the exclusion condition of states at ^ 1. Therefore, we 
shall numerically clarify the phase diagram in the two- 
dimensional model-parameter space. 

Let us consider A with M rows of L sites wrapped 
on a cylinder and take an even number L and the limit 
M 00. Then, the site j e A is specified by I E 
and m € [1, M]. For this system, we can define the trans- 
fer matrix T(L) connecting the NNN pair of rows, and 
denote its eigenvalues as {Xa{L)} (a specifies a level). 
The conformal field theory provides direct expressions 
for c and Xa (the scaling dimensions) of critical systems 
by using the eigenvalues 0, : 

-Cln[A/(i)] ~L/-7^c/6L + 6/L^ (3) 
~On[Xa{L)]+ Cln[\i{L)] ~2nxa/L, (4) 

where A/ is the largest one corresponding to the ground 
state. C,, /, and b are the geometric factor (i), a free 
energy per site, and a nonuniversal constant, respec- 
tively. Here, it should be noted that the discrete sym- 
metries of the lattice Hamiltonian, e.g., the translation 
T (e.g., ai^rn m+2,m); the space inversion V (e.g., 
CT/,m (JL-i+2,m), and the S3 symmetry [cr;,m 9{(^i,m) 
for g S S3] are crucial not only for a reduction of compu- 
tational efforts, but also for the proper specification of the 
level a This can be clearly demonstrated by evalu- 
ating the third scaling dimension xp &i {K, V) = (00, 0). 
According to Ref. 0], the staggered polarization opera- 
tor is invariant for all g € S3 and has the wave number 
7r/a so that the corresponding level — Cln[Ap(L)] should 
be found in the subspace specified by these symmetries. 
In Table we give the numerical data for the scaling 



TABLE I: The size dependence of the scaling dimension xp. 
The fitting xp{L) — xp(oo) + bi/L^ + h2/L* is performed 
using the data of L =10, 12, and 14. 



L = 8 


10 


12 


14 00 


xp{L) 1.6801468 


1.6065435 


1.5710939 


1.5510358 1.5007 



dimension xp{L) estimated from Eq. Q). The extrap- 
olated value does not deviate more than 0.1% from the 
theoretical one, |, and the level is not the third one in 
the whole space, so the characterization of the excitation 
levels is essentially important (for other exponents, see 
Ref. [13). 

Now, according to the prediction of the crossover, 
it is plausible that the critical RG flow starts out of 
[K, V) = (00, 0) and arrives at the point (0, V^^). Thus, 
at this stage, our numerical task is to determine the 
line V3(A'). For this purpose, we shall employ the phe- 
nomenological RG (PRG) method Let us denote the 
left-hand side of Eq. as /S.Ea{K^ V, L) (i.e., an excita- 
tion gap), then we shall numerically solve the following 
PRG equation for a given value of K with respect to V: 
LAEa,{K,V,L) = L'AEc,{K,V,L'). Since this is satis- 
fied by the gap AEa{K, V,L) cx the obtained value 
can be regarded as the size-dependent estimate of the 
transition point, say V3{K,L) [we take L = {L + L')/2 
and i' = L — 2 in the following]. Further, there arc two 
critical fixed points connected by the RG flow, so a rela- 
tionship between lower-energy excitations on these fixed 
points, namely the ultraviolet-infrared (UV-IR) operator 
correspondence, is quite important for the choice of the 
excitation a. Along the flow, the conjecture, Sj —^ Sj 
and S* s* , has been proposed ^3 • Sj is the magnetic 
operator on the IR fixed point with the three-state Potts 
criticality, whose scaling dimension is xg ^ j^. Since the 
excitation Sj provides the lowest energy level, we shall 
focus our attention on A£'j(Ar, V, L) stemming from the 
level of Sj on the UV Gaussian fixed point. 

The exact diagonalization calculations of T(L) with 
L =4-14 are performed by using the Lanczos algorithm, 
where the discrete symmetries T, V, and g e S3 are uti- 
lized. We plot examples of L and V dependences of the 
scaled gap {L/2t:)AEs{K,V, L) in Fig. HI and find the 
crossing points. While non-trivial finite-size corrections 
may affect their behaviors, we shall extrapolate them to 
the thermodynamic limit according to the finite-size scal- 
ing argument [i^: Suppose a single-power formula, i.e., 
V3{K, L) — V-i{K) (X L~'^^ , then the exponent is given as 
-03 = C1J3 + l/i/, where lo^ and v are the correction expo- 
nent and the critical exponent of the correlation length. 
For the F three-state Potts model, a;3 = | Q and 
1/1/ = |, so we shall use the formula with ijj^ = 2 and 
extrapolate V3{K,L) to the limit. 

For convenience, we define the reduced couplings 
(m, u) = (e~^,l — e~^) and compactify the parame- 
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FIG. 1: The V dependence of {L/2-k)AEs[K,V, L). From 
left to right, K = 2.0, 1.0 and 0.5, respectively. The cor- 
respondence between marks and system sizes is given in the 
figure. The crossing points give V3{K,L). 




FIG. 2: The phase diagram of the Hamiltonian Q in the 
space of the reduced parameters {u,v) — {e~^,l — e^^). 
Open circles (double circles) with a curve exhibit the phase 
boundary 113 (w) [t;2(u)], which belongs to the three-state Potts 
(Ising) universality class (arrows show the directions of the 
RG-flows). The filled triangle on the v biaxis denotes the F 
three-state Potts critical point (1,113''), the filled (half- 
filled) circle corresponds to the ground state of the AF three- 
state Potts model on A (A_) with the Gaussian criticality. 



ter space within a unit square region. Then we draw 
the phase boundary line v^^u) corresponding to V3{K) 
(open circles with a fitting curve) in Fig. [21 The size- 
dependent data V3{K,L) are also given by other marks. 
From this figure, we can find the following: The ex- 
trapolated boundary line starts out of the filled circle 
(u, v) = (0, 0) linearly with the increase of u, which agrees 
with the crossover argument, i.e., V3{K) cx e"^"^^ with 
(j)3 = {2-xp)/{2-x^) = 1 Hiil. Then, V3{u) monoton- 
ically increases and finally terminates at the filled trian- 



gle {u,v) = {l,vf^) with vf^ := 1 — e~^3 , as expected. 
Consequently, the boundary line connects the Gaussian 
and the three-state Potts critical fixed points, which is in 
accord with Zamolodchikov's c-theorem 0|. 

The line V3{K) is now known to separate the disor- 
dered and the partially ordered phases, because the spin 
degrees on A_ are still disordered, at least, on the decou- 
pling line, K — and V > V^^. Therefore, another phase 
boundary to the ordered phase exists. To see this, let us 
investigate our model for large enough V. Since the F 
order is well established on A+, we can approximately re- 
place the spin degrees on A_|_ by a spontaneously favored 
value, e.g., r. Then, the subsystem on A_ is effectively 
described by the AF three-state Potts model under the 
uniform magnetic field; the effective Hamiltonian is 



(5) 



(j, fc e A_). Again, the AF Potts criticality at [K, V) = 
(0, 00) is perturbed by two relevant operators, but unlike 
in the case of Eq. (j2J) the S3 symmetry is absent, so that 
a different type of crossover may be expected. For large 
K, since aj cannot take r, the subsystem is effectively 
described by the Ising model. Racz and Vicsek investi- 
gated the Hamiltonian ^ by the use of the MC method, 
and obtained the boundary line |2^, so that our next 
numerical task is to determine the second line V2{K) to 
complete the phase diagram of our Hamiltonian 

The relevant level in our PRG calculations may corre- 
spond to the magnetic-type excitation on the Ising criti- 
cal fixed point AEcr{K, V, L) with the scaling dimension 
Xg. = i, and it can be found in the above of the levels 
which are to be degenerate to the ground state. So, sim- 
ilarly to V^lK.L), we can calculate the size-dependent 
estimates of transition points V2 {K, L) from the cross- 
ings of the scaled gaps LAE„ {K, V, L) . Then, extrap- 
olating them to the thermodynamic limit by u sing the 
formula V2iK,L) - V2{K) cx with V2 = 3 HI, we 

obtain V2 (K) [the corresponding value V2 (u) is exhibited 
by the double circles with a curve in Fig. |2] . Although 
Eq. jSJ becomes exact in the limit V ^ 00, the obtained 
boundary line is qualitatively comparable to that in Ref. 
|26| also for finite couplings. Here, let us examine the 
limit V2{K 00). Suppose that the complete F or- 
der is established on A_|_, then the limit is given by the 
transition point of the Ising model on A_, V^^ ~ 0.8814 
(1 — e~^2"' ~ 0.5858 in Fig. EJ. However, the spin con- 
figurations which are neglected in this argument exist in 
the original Hamiltonian ((SJ, so V2(oo) may deviate from 

^2(00) ~ 0.8820, is quite 
This implies that the contributions from 
the neglected spin configurations may be small around 
¥2^, and thus Eq. © is expected to provide a very good 
description even near V2^. 

Now, in order to evaluate the conformal anomaly num- 
bers, we perform three point fitting on Eq. Q based on 



¥2^. Our numerical estimate 
close to 
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TABLE II: The size-dependent estimates of the conformal 
anomaly numbers Cq{L) at {K, V) = (1, Vq(l)) (g =2, 3). The 
fitting Cq{L) — Cq{oo) + fci/L^ is performed using the data of 
L =10 and 12. 



L = 6 


8 


10 


12 


<x 


C2{L) 0.5134201 


0.4359578 


0.4399469 


0.4611012 


0.5092 


C3{L) 0.8012989 


0.7682470 


0.7714379 


0.7815416 


0.8045 



L — 2, L, and L + 2, and estimate the size-dependent val- 
ues Cq{L) (q =2 and 3). As an example, we show them 
at intermediate coupHngs {K,V) = {l,Vq{l)) in Tabic 
im Although nonmonotonic size dependences are visi- 
ble, the extrapolated data agree well with the theoretical 
values, ^ and |. For parameters close to the IR fixed 
points, i.e., (iC, V2iK)) with K > 1 and (A", VsiK)) with 
K <^ 1^ there are only small L dependences, and it is 
easy to confirm the criticalities. However at {K,Vq{K)) 
far away from the IR fixed points their precise estima- 
tions become difficult for the present system sizes. 

Finally, we make a remark on the behavior around 
(K, V) = (0, oo) [i.e., the half-filled circle (u, v) = (1, 1) in 
Fig. 12 . We have considered the model (0) possessing the 
global S3 symmetry. With the increase of V, S3 is broken 
at V3{K) and reduced to the Z2 symmetry with respect 
to the interchange of unfavored two spin states ( "charge 
conjugation") for V > Vs{K). While the lattice Hamil- 
tonian (jSJ well describes the transition accompanied by 
the Z2 symmetry breaking |26l l27l|. corresponding field 
theory in the scaling limit is required to clarify, for in- 
stance, the phase boundary line around {K, V) ~ (0, 00). 
For this issue, it is plausible that the symmetry-breaking 
negative field (oc K) couples with the uniform magneti- 



zation Sj and brings about a competition to the energy 
operator ej (a coupling oc e~^). By borrowing the sine- 
Gordon expression in Rcf. Q, this competition may be 
formulated by the action preserving the sublattice sym- 
metry and the charge conjugation but breaking the cyclic 
permutation symmetry as 



d^a; 



^ cos Pip + 11 cos ( '^V 



(6) 



where Ac=i is the free boson action, and p is the dual 
field of p. The possibility that the dual sine-Gordon the- 
ory © describes the Ising criticality was mentione d by 
Delfino in relation to the Ashkin- Teller model [11 lisf . 
Further, note that the action © can become self-dual at 
0^ = 47r which is known to exactly show the Ising crit- 
icality 01 . Based on this expression, we can predict a 
shape of the boundary line as follows: K oc e~^^^^'>'^^ 
with 02 = (2 — Xs)/{2 — x^) — |. This seems to coin- 
cide with the rapid change of the boundary observed in 
Fig. [5] However, the precise estimation of 02 is outside 
of the scope of our present research; the investigation to 
confirm this prediction is now in progress. 



To summarize, we have investigated two crossovers of 
criticalities in the square-lattice three-state Potts antifer- 
romagnet with staggered polarization field. On the basis 
of the dual sine-Gordon field theories, we have also given 
an argument on the criticalities and the shapes of the 
phase boundaries. For more detailed study, we are now 
performing the Monte Carlo simulation calculations to 
show the evidences for supporting the present result; we 
will give them in the forthcoming article. 

Main computations were performed using the facilities 
of Yukawa Institute for Theoretical Physics. 
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